A simple approach to the calculation of retarded dispersion forces 
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Abstract 

We propose a phenomenological Hamiltonian for the interaction of neutral macroscopic bodies 
with the electromagnetic field. Subsequently we revisit the assumption according to which the 
retarded interactions between neutral macroscopic bodies can be obtained through an additive 
principle i.e. summing the volume elements defined by the macroscopic bodies. 
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I. INTRODUCTION 



As it is very well known the discovery and the first calculation of the retarded forces 
between macroscopic bodies, due to the fluctuations of the electromagnetic field, can be 
found in Casimir's 1948 seminal paper A few years later the Lifshitz theory appeared. 
Within the frame of this theory a powerful expression describing the force at finite tem- 
perature between two semiinfinite dispersive media separated by a slab was obtained j^. 
However the fundamental method of calculating these interactions at zero temperature is 
still due to Casimir. One has to single-out the electromagnetic field normal modes con- 
sistent with the boundary conditions imposed by the presence of the macroscopic bodies. 
The shift of the zero point energy of the field with respect to vacuum is then considered as 
the sought interaction energy. In the case of complicated geometry the evaluation of the 
Casimir force is very difficult, therefore various approximation methods were developed such 

n nn 

as the so called proximity-force approximation (PFA)[3|], the semiclassicalj^, |5j and optical 
approximatiojisj^] and path-integral methods [3]. We are interested here in an approach used 



in reference |8| and discussed by Mostepanenko and Sokolov[9[]. The interested reader can 

lOl. It is assumed 



find a review of the subject in the comprehensive report by Bordag et al. 
that the dependence of the interaction potential of macroscopic bodies on distance can be 
determined by summing the interaction potential between single atoms (molecules) belong- 
ing to different bodies. The coefficient of the expression so obtained, which depends on the 
kind of material involved, is renormalized in order to be in agreement with the interaction 
between flat plates of the same material as derived by the Lifshitz theory. The afore men- 
tioned additivity principle works very well in the case of low density media though it is not 
said that the principle could be extended in general and particularly in the case of metals. 
However because of the simplicity and usefullness of this approach we are here revisiting 
it with the idea according to which the retarded forces between two macroscopic bodies 
might be understood as being due to the attraction in pairs between their volume elements. 
Evidently such volume elements must be considered small as far as the differential calculus 
is concerned, yet rather large compared to interatomic distances. We will delineate our idea 
proposing a phenomenological Hamiltonian for the interaction of the electromagnetic fleld 
with a macroscopic body whose properties are phenomenologically described. 

In the following sections we will discuss the proposed Hamiltonian showing where it 
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leads. Subsequently we will apply our results to a few geometrical dispositions of simple 
macroscopic objects. 



II. A PHENOMENOLOGICA PERTURBATION HAMILTONIAN 

Consider two neutral, dielectric and homogeneous bodies in fixed positions. The presence 
of such bodies will modify the Hamiltonian of the electromagnetic field which will be written: 

n = no + n' (i) 

where Ho is the Hamiltonian of the free electromagnetic field and where H' indicates the 
perturbation due to the presence of the two bodies. Assuming as unperturbed field state the 
vacuum we will calculate the energy shift due to the extra term of the Hamiltonian. The part 
of the shift which depends on the configuration of the two bodies will be interpreted as their 
potential energy at zero temperature. We consider and propose now a phenomenological 
perturbation as follows: 

n' ^ -I- I P -EdV -I I M-UdV (2) 

2 Jv 2 Jv 

where E, H are the electric and magnetic fields, P and M are the polarization and mag- 
netization vectors. The integral is extended to the volume occupied by the bodies. Always 
within a phenomenological approach the Fourier components Pk and Mk may be considered 
connected to those of the electric field and magnetic field by a linear relation. 
We set: 



P = iV2^J2^(^(^) (ok.Ae^'"'' - 4,Ae"'^"") eA(k) 
M 



k,A 

= iV2^ J] \/^7(^) («k,Ae^''-^ - <Ae-^'' '") ^ ^^^^^ 



(3) 
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k,A 

where ak, Ok annihilation and creation operators and eA(k) indicates the polarization. 

For the field we have used only transversal waves and we set to unity the side of the cube 
which defines the periodic conditions. The P{k) and 7(A;) coefficients are defined as: 

3 e{k) - 1 



m 



AT^e{k) + 2 
3 ^i{k) - 1 

47r/i(A;) + 2 



(4) 



where e and jj, are the electric and magnetic permittivities which expresses the usual 
Clausius-Mossotti formula. We then obtain: 



eA(k) • e^(k )[3{k) + -^[k) 



k k' i 

a\L,\a\i' ,vVi{\<i + k') - aj,_Aakv^^i(k' - k) - ak,Aaj,,^t;i(k - k') + 4,;^aj^,^fi(-k - k') (5) 
The index i = 1,2 indicates the summation over the two bodies. We furthermore defined: 

^kiv) = [ e'^^-^'dV, (6) 



III. THE INTERACTION ENERGY AT LARGE SEPARATIONS 

We stress that we consider the perturbation Ti' above to hold only for the calculation of the 
interaction energy of neutral macroscopic bodies at large separations. Let d be the shortest 
distance between the surfaces of the two bodies. The expression "large separation" usually 
means d » Xi where Aj are the absorption wavelengths of the material under consideration. 
The mutual energy of the two bodies at zero temperature is obtained by perturbation at 
first order. In the case of large separations the calculation simplifies greatly because 
and 7(/c) can be safely substituted by /5(0) and 7(0). In fact the major contribution to 
the integral below comes from the k < 1/d region where d, as previously mentioned, is the 



shortest distance between the surfaces of the bodies 
Thus we asymptotically obtain: 

fic 

Ui,2 = m0)f32{0) + 71(0)72(0)] X 
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Setting q = k + k' the expression ([7]) becomes: 



He 

Ui,2 = mo)m) + 71(0)72(0)] X 



/ .c(.(c)..(-q) . e.c.) / .k^^ (1 . ii,^) (8) 



The integral over k because of a lack of a quenching factor obviously diverges. Yet within 
the scope of our calculation we can neglect the infinities and consider only the singular part 
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in the origin (g = 0) which dominates the asimptotic behaviour. The present procedure is 
analogous to the usual renormalization techniques which make use of a cut-off. In our case 
the divergent part of the integral is already separated without ambiguities. Furthermore this 
part can be ignored because it does not depend on the system configuration and therefore 
it can not give origin to forces among the considered bodies. The final result, indicated in 



expression ([9]), can be obtained rather laboriously following the indications of reference |12| : 

[/i,2 = /icy4i,2 j dc[{vi{q)v2{-q) + c.c.)q^\nq (9) 

where: 

23 

Ai,2 = mo)m) + 71(0)72(0)] (10) 

For two perfect metals [^(O) = 00, /i(0) = 0] the constant Am,m is: 

_ 1035 1 

Am,m - 3g^Q ^2vr)5 ^ ^^^^ 

Now using the tables of the Fourier transforms of the generalized functions [131 we have: 

rfqgMnge^^-^ = -^^^ (12) 
7rr' 



hence recalling ([6]) we obtain: 



Ui,2 = -^c5i,2 / c/ri / dr2, ^ (13) 

JVi JV2 I'^l ~ ^2 I 

where: 

23 

Bi,2 = imw) + 71(0)72(0)] (14) 

In the case of two perfect metals the constant is: 

_ 1035 1 

-Dm,m — -^737 n (loj 

256 TT'^ 

As it can be seen from expression f|T3|) our results confirm Mostepanenko's assumptions!^ 
yielding in addition the coefficient as a function of the characteristics of the material. 



IV. A FEW EXAMPLES 

In the case of two alike homogeneous dielectrics of very low molecular density and where 

the molecules are assumed of the same kind 7(0) = and /5(0) becomes Na{Q) where is 

the molecular density and a(0) the molecular polarization. From (1131) we then have: 

23 r C 1 

U = hcai(0)a2(0) / Ndr^ / Ndr2-, ^ (16) 

47r J J |ri-r2|^ 



Which means that the interaction potential of the retarded Van der Waals forces between 
two molecules is: 

09. 1 

(17) 



which is a well known result. 
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23 1 

— /icai(0)Q;2(0)i 7z 

4,n |ri-r2|^ 



li 



Shall we move on now to the classic case of two parallelepipeds made of perfect metal 
whose volumes are so defined: 



L L 

< Xi < - 

L L 

<X2< - 

2-2 



The volume Fourier transforms are: 



L L 

-<yi< - 

L L 

— _ ^2 _ — 



wi(q)^2(-q) = 16 



sm 



qyL 



sin ^ \ (1 



- D < ^1 < 
a < Z2 < a + D 



(18) 
(19) 



1 



(20) 



To integrate (!2Ql) we find it convenient, not to mix up the calculations, to write: 



sm 



sm ^ sm 



Qx I Qx Qx 

regarding the first of the factors on the right (assuming L large) as 

qxL 



(21) 



sm 



Qx 



T^K^lx) 



(22) 



we then calculate the limit U ^ L once the integration has been carried out. We follow the 
same route with the factor depending on qy. We are then led to the integral: 



(23) 



which can be evaluated using again reference 13|. One then finds: 



2(27r)^L^ - + 



1 



(24) 



{a + Df ia + D)\ 

Inserting the previous result into expression we obtain the interaction energy for two 
metallic parallelepipeds at zero temperature: 

1 2 



1035 fc^2 / 1 
7680 7r2 



+ 



[a + 2Dy (a + D)3 



(25) 
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which at the hmit D ^ a becomes: 



1035 hcL'^ HcL^ 



-0.01365^^ (26) 



7680 a3 ■ a3 

Expression (!26l) has to be compared with Casimir's more elegant resuh [l|: 



vr^ HcL'^ _ hcL 
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0.01370^^ (27) 



720 a3 ■ a3 

The small discrepancy between the two results deserves a brief comment. Though its origin 
is not clearly evident, we think that it comes from the fact that our Hamiltonian, referring 
to the bulk of the surrounding bodies, does not produce an exact description of the field in 
the proximity of the surface where the Clausiu-Mossotti formula may be inadeguate. In the 
Casimir's approach the appropriate boundary conditions are generally referred tho the field, 
thus the description of the situation on the surface is certaninly more exact. 

Consider now the interaction u of a metallic corpuscle of infinitesimal volume dV at 
distance y from a metallic hemispace. From expression (fT3l) with B12 = Bm,m and using 
cylindrical coordinates we have: 

/■27r roo roo 

u = -HcBm,mdV / d(f pdp / dz (28) 

Jo Jo Jo [p^ + {z + y)"^] ' 

which is: 

u = -hcBm,m^^dV (29) 

Using this last expression one can calculate the interaction between a very thick and infinitely 
extended metallic plate and a generic body. Let this last be a metallic sphere of radius R 
at a distance r = d + R. We obtain as interaction energy: 

^ /•2n pI pR P^dp 

v,s = -ncB.^,^-j^ d^jjcosej^ (7T7^ 

( R _^ 1 R \ 



with: 



Bm,n.— = 0.04289 (32) 



This result should be compared with that of Balian and Duplantier 
tering approximation for d << R is: 



15| which in two scat- 



where: 




0.03978 



•) 



We think that, apart from a few shght differences in the numerical coefficients, our expression 
can lead to a good approximation with respect to the exact calculation. 

Although the dependence of the Casimir effect on the temperature is not within the scope 
of the present brief note, we think appropriate to mention at least the consequences that 
our Hamiltonian implies on the subject. Our Hamiltonian includes diagonal terms which, 
for each body which may be present, generate a field energy equivalent to the body volume 
times the black body radiation energy density as shown by Stefan-Boltzmann's expression. 
But this is still a well known result. It is more interesting to investigate on the effect of 
temperature on the dispersion forces. To this end it is sufficient to insert as a factor of 
the integrand in equation ([7]) the quantities n{k\ T) + n{k, T) where n(A;, T) is the average 
number of photons of wave number k at temperature T obtained from Planck distribution. 
The expression which comes out is not analytically menageable, however the insertion of 
the above mentioned factor introduces a natural cut-off which implies that the quantity 
Ui,2{T) can be easily evaluated using numerical methods. However an analytical result can 
be obtained in the case of high temperatures developing to the first order n{k, T) and setting 
it equal to KT/{fick). This procedure introduces a k factor at the denominator of equation 
([7]). Then from a power counting one can see that in this limit the interaction f/i 2(r) 
between two volume elements at distance r will behave as As a consequence in the 

case of two perfect metals in Casimir's configuration we will obtain, instead of expression 

nn 

( l26ll . an interaction energy proportional to —KT{L'^ /a?') which is a well known result |l6Ul7l|. 
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